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ABSTRACT 

We studied the structure of two old open clusters, NGC 1245 and 
NGC 2506, from a wide and deep VI photometry data acquired using the 
CFH12K CCD camera at CFHT. We devised a new method for assign- 
ing cluster membership probability to individual stars using both spatial 
positions and positions in the colour-magnitude diagram. From analyses 
of the luminosity functions at several cluster-centric radii and the radial 
surface density profiles derived from stars with different luminosity ranges, 
we found that the two clusters are dynamically relaxed to drive significant 
mass segregation and evaporation of some fraction of low-mass stars. There 
seems to be a signature of tidal tail in NGC 1245 but the signal is too low 
to be confirmed. 
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1 INTRODUCTION 

Open clusters are useful examples for studying the dy- 
namical evolution of stellar systems and the chemical 
evolution of the Galactic disc, because the stars in 
an open cluster are assumed to have the same chemi- 
cal and dynamical properties owing to of their shared 
origin. Because open clusters orbit about the centre 
of the Galaxy near the Galactic plane, they are fre- 
quently disturbed by both giant molecular clouds and 
the tidal forces originating from the Galactic disc. In 
particular, old open clusters can be tracers of the 
structural and evolutionary history of the Galactic 
disc (Friel 1995). 

The structure of the old open clusters can be 
changed by dynamical evolution driven by internal 
and external tidal forces. Energy equipartition due to 
encounters among member stars (lSpitzeiill940l ) leads 
to mass segregation, which is characterized by the con- 
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centration of high-mass stars in the cluster's centre, 
and an isotropic distribution of low-mass stars ram- 
bling throughout the cluster. Mass segregation results 
in the steeper slope of the mass funct ion in the outer 
region of the cluster (jAnn fc Leell2002l ) . Some fraction 
of low-mass stars evaporates from the cluster to create 
the halo surrounding the cluster ([Eggen 1993). The 
evaporation of low-mass stars also increases the size 
of the cluster because of the reduction of the binding 
energy. 

Although internal interactions among member 
stars make the distribution of stars more spher- 
ical, perturbation by Galactic tidal forces results 
in an elongated shape. The elongated morphology 
of open clu sters has been reported for the Hyade s 
(|Oortlll979l '). Pleiades ("Raboud fc M ermilliodlll998al'). 
NGC 2287, and NG C 2548 (Bcrgond. Teon fc Guibera 
I2OOII ). IChen. Chen and Shul (|2004l ) also reported the 
elongated structure of 31 open clusters, which includes 
young star clusters that are a few million years old. 
ISharma et all ([2009) reported the elongated corona of 
6 open clusters among 9 open clusters. They could not 



© 2013 RAS 



2 Lee, Kang & Ann 



Table 1. Physical parameters of NGC 1245 and NGC 2506 taken from Paper I 



Cluster 
name 


a 
J2000.0 


5 
J2000.0 


I 
(dog.) 


b 

(deg.) 


E(B - V) 


(y-A/„)o 


age(Gyr) 


[Fe/H] 


NGC 1245 
NGC 2506 


03 14 47 
08 00 02 


+47 13 53 
-10 46 15 


146.65 
230.56 


-8.93 
+9.94 


0.24 
0.035 


12.25 
12.47 


1.08 
2.31 


-0.08 
-0.24 



Notes. - Units of right ascension are hours, minutes, and seconds; units of declination are degrees, arcminutes, 
and arcseconds. 



find a correlation between the age of a cluster and its 
core shape. 

One possible origin of the elongated morphol- 
ogy of open clusters is a stretching of the stellar 
distr ib ution along the cluster orbit (|Chuinak et al.l 
bold ) . iDavenport fc Sandauisti (|2O10l ) represented the 
halo morphology of M67 is elongated shape us- 
ing SDSS data, which is roughly aligned with the 
proper motion of the cluster. A relevant structure 
is the tidal tail caused by the Galactic tidal field 
(|Chumak fc RastorguevI 12006) . Because the evapora- 
tion of low-mass stars is likely to be isotropic, if there 
is no external force, the elongated shape of the halo of 
open clusters seems to be caused by the same mecha- 
nism that creates the tidal tails. Thus, it seems likely 
that a detailed analysis of halo structure provides im- 
portant information about the dynamical history of 
open clusters. However, analysis of the spatial struc- 
ture of open clusters is not easy, because it is ex- 
tremely difficult to assign cluster membership without 
having kinematic information. 

The membership of cluster stars is critical for 
any investigation of the dynamical structure of clus- 
ters. Earlier studies of dynamical evolution and halo 
structure were carried out for nearby clusters whose 
member stars were well defined by proper motion 
studies ('Mermilliod et al."l990|: |Perrvman et al.ll 19981 : 
iRaboud fc Mernulliod 1998a ;lj). Since the kinematic 
data for membership determination is difficult to be 
attained for most of open clusters, especially for the 
faint stars, the colour-magnitude diagram (CMD) of 
clusters and surrounding field regions are used to cor- 
rect the field star contamination in th e cluster lu- 
minosity function (|eg.. Sung et al.lll999l ). There have 
been some attempts to determine membership proba- 
bility of individual stars statistical ly using spatial den- 
sities (IChen. Chen and Shul 20041 ) or stellar densities 
in CMD i Janes fc HoaluOllI ). Here, we introduce a 
new method for assigning cluster membership proba- 
bility based on locations in the CMD and projected 
spatial positions in the sky. 

The aim of the present study is to understand 
the structures and dynamical properties of two old 
open clusters, NGC 1245 and NGC 2506, for which 
VI CCD photometry data has been reported by 
iLee. Kang fc Ann! ()2012h (Paper I). Because the pho- 
tometry was deep and wide enough to analyze the 
surface density distribution of low-mass stars around 
the tidal radfi of NGC 1245 and NGC 2506, we expect 



to be able to describe the halo structures of the two 
clusters. We also expect to see the effects of dynam- 
ical evolution on the structure of open clusters when 
analyzing the radial dependence of LPs. 

In §2, we describe the observational data and the 
cluster parameters. We describe membership probabil- 
ity in §3. The properties of dynamical evolution and 
halo are analyzed in §4. A brief summary and conclu- 
sions are given in the final section. 



2 OBSERVATIONAL DATA 

The data used in the present study are from the deep 
and wide field VI CCD photometry observations we 
conducted in a previous study (Paper I). A detailed 
description of the observations and data reduction are 
given in Paper I; however, we briefly describe the pho- 
tometry here. We obtained the data using a CFH12K 
mosaic CCD camera mounted at the prime focus of 
the CFHT. The CFH12K CCD camera has the field 
of view 42' x 28' with a pixel size of 0."206. We ob- 
served 6 regions in NGC 1245 and 3 regions in NGC 
2506 (all of which were centred on the cluster centres), 
and which covered 84' x 82' and 42' x 81' for NGC 1245 
and NGC 2506, respectively. 

We followed the standard procedures for CCD re- 
duction using IRAF and carried ou t PSF photome- 
try using DAOPHOT JStetsonlllOSTJ ). The seeing size, 
defined as the FWHM of the stellar images, varied 
between l.'l and 1.'6 from region to region. For cal- 
ibrati ng our photom etry, we observed the SA 98 re- 
gion (|Landoltlll992l ) on the same night. The limiting 
magnitude of the photometry is V ~ 23. The cluster 
parameters are listed in Table 1; they were derived 
using the CMDs as described in Paper I. 



3 MEMBERSHIP PROBABILITY 

Traditionally, proper motion study is the most reli- 
able method for determining the membership criteria 
of cluster stars. However, this method requires two 
observations separated by a long time interval. There- 
fore, most previous studies utilized membership crite- 
ria based on the locations of cluster stars in the CMDs 
(|eg.. Sung et al.lll999 ). This approach is very simple 
and effective when the main sequences do not overlap 
with wedge patterns of field stars in the CMD. 

The main sequence and giant stars of a cluster are 
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Figure 1. Color-magnitude diagrams of NGC 1245 and NGC 2506 with the stars in the central region {R < 10'). Green 
solid lines in the panels indicate the photometric membership boundaries of the cluster stars. The blue dashed line indicates 
the rough positions of the main sequence to divide the non-member stars into those above the main sequence and those 
below the main sequence. 
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Figure 2. The normalized spatial density distribution of field stars which arc distributed in the outer region of the 
membership boundary in Fig. 1. Upper pannels are for NGC 1245 and lower pannels are for NGC 2506. The left pannels 
and right pannels show the distributions of the field stars located above and below the blue dashed lines in Fig. 1, respectively. 
The middle pannels represent the density distribution of all field stars which are located outside the membership boundaries 
in Fig 1. The colors represent the density increment of Ag(x,y) = 0.14 from the lowest level as blue color and the highest 
level as red color with three intermediate levels, cyan, green and yellow among which green color represents the density 
level of g{x, y)=l. 
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Figu re 3. Correlation betwe en proper motion member- 
siiip l lCtiiu fc van AItenalll98ll ) and photometric statistics 
for stars in NGC 2506. Stars between dotted lines have 
membership differences smaller than 20%. 



distributed along the cluster isochrone in the CMDs, 
whereas non-member stars are primarily distributed 
in a wedge shape with a random distribution outside 
the wedge pattern. There is also a clear difference in 
the spatial distributions of the cluster member stars 
and those of non-member stars. The cluster stars are 
likely to be centrally concentrated, whereas field stars 
are uniformly distributed, which leads to an outward 
decrease of the ratio of the number of member stars 
to the number of field stars. Thus, it seems plausible 
that we can differentiate the cluster members from the 
field stars using CMDs and the spatial distributions of 
stars in the cluster region. 

According to the stellar evolution theory, clus- 
ter members are distributed along the main sequence 
and giant branches in the CMD. To determine the 
boundaries of the two clusters within which the mem- 
ber stars are supposed to be located in the CMD, we 
derived the density distributions across the main se- 
quence at each magnitude bin using stars in the cen- 
tral region, R < 10' and fitted Gaussian functions to 
the density cross-sections along the main sequence. 
We smoothly connected the locations correspond to 
the 3 a density levels to determine the photometric 
membership boundary in the CMDs (Fig. 1). For the 
giant stars we defined the boundaries based on the 
visual inspection of the giant stars in the CMD. We 
neglected the blue stragglers in defining the member- 
ship boundaries. The magnitude limits are 1/ ~ 23 for 
the clusters. 

We define the membership probability of a star 
i, Pi, as a function of the coordinates in the four- 
dimensional (4-D) space that consists of two spatial 



coordinates in the sky and two photometric parame- 
ters constituting the CMD, namely x, y, V, V—I. Thus, 
we consider a 4-D membership probability as follows. 
We counted the number of stars in a unit volume 
of the 4-D space centered on a point {x, y,V,V — I) 
and then corrected for the contribution by field stars. 
We derived the field star distribution in CMD, i.e., 
Nf{V, V — I), to correct for field star contamination 
using the stars in the outermost regions R > 40' for 
NGC 1245 and R > 35' for NGC 2506. If we assume 
that the distribution of field stars is homogeneous in 
the cluster region, the distribution function of cluster 
member stars in the 4-D space, Nc{x,y, V,V ~ I), is 
determined from the following equation: 

N,{x,y,V,V-I) = N{x,y,V,V-I)~Nf{V,V^I){l) 

where N{x, y,V,V ~ I) is the distribution function of 
stars in the 4-D space, including both member and 
non-member stars. 

Because field stars are mostly disc stars which 
are located in the Galactic plane where considerable 
density variation is expected, it is better to take into 
account the density variation of field stars. To inves- 
tigate the variation of the spatial distribution of non- 
member stars, we derived the spatial distribution of 
field stars, using stars outside the cluster boundary, 
by calculating the normalized field star density distri- 
bution g{x,y) defined as 



g{x,y) = n{x,y)/ < n{x,y) > 



(2) 



where n{x, y) is the spatial number density of stars 
outside the photometric membership boundaries in 
the CMD and < n{x, y) > is the average value of 
n(x,y). Thus, the mean value of g(x,y) is equal to 1. 
We derived n{x, y) by counting the number of stars 
in the smoothing box of 4' x 4', with sampling in- 
terval of 2' X 2'. These sampling size and sampling 
interval were used through out the paper for the 
derivation of the surface number densities such as N 
(x,y,V,V-I) and g(x,y). The sampling size and sam- 
pling interval in the CMD are AV' = 1.04,0.52 and 
/\{y-I) = 0.269, 0.134. We used a Gaussian smooth- 
ing algorithm in 4-D space to derive the surface num- 
ber densities. 

Fig. 2 shows the spatial distribution of non- 
member stars, represented by g{x,y). The distribu- 
tions of field stars below and above the main sequence 
are plotted in the right panels and left panels, respec- 
tively. The middle panels show the spatial distribution 
of all the field stars in the CMDs. The upper panels 
are for NGC 1245 and the lower panels are for NGC 
2506. There are some differences in the spatial distri- 
butions of the stars in the right panels and left panels. 
In particular, the stars in the upper right panel, i.e, 
field stars below the main sequence of NGC 1245 show 
higher density in the middle south region of the clus- 
ter while those above the main sequence of NGC 1245 
show higher density in the north-east part of the clus- 
ter. We suppose that this difference is caused by the 
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Figure 4. Spatial distribution and CMD of an artifi- 
cial cluster generated by Monte-Carlo simulation for NGC 
1245. 

gradient in the density distribution of field stars to- 
ward NGC 1245. The stars above the main sequence 
are likely to be foreground stars while those below the 
main sequence are likely to be background stars be- 
cause the stars above the main sequence are brighter 
than those below the main sequence if they have sim- 
ilar colors. 

Using g{x,y), we rewrote the 4-D distribution 
function of cluster member stars Nc{x, y,V,V — I) as 

N,ix,y,V,V-I) = N(x,y,V,V^I) 

^ Nf{V,V ~ I)g{x,y) (3) 

Then, the membership probability, Pi, of a star i is 
defined as 



We derived the membership probability of all stars in 
the observed regions of the two clusters, however, we 
assigned Pi = to the stars outside the photometric 
membership boundary in the CMDs (Fig. 1). 

To check the validity of the derived values of Pi, 
we compared the membership probabilities of stars 
in NGC 2506 derived from the present me thod using 
P{x, y, V V - I) and P„ of (IChiu fc van Altena 1981f ) 
in Fig. 3. IChiu fc van Altenal ( 19811 ) derived the mem- 
bership probabilities P^ of 724 stars in NGC 2506 
using two photographic plates with an epoch differ- 
ence of 56.9 years. Their sky coverage is only 5. '3 in 
radius, with a magnitude limit oi V ~ 16. As can 
be seen in Fig. 3, the majority of stars with a high 
membership probability (P,, > 0.7) also have a high 
P{x, y,V,V — I) value. However, a significant fraction 
of high P{x, y,V,V — I) stars have a P^ less than 
0.7. This is because the proper motion study was con- 
fined to the stars in the central regions of the cluster, 
where high values of P{x, y,V,V — I) are assigned to 
stars if their CMD positions are within the photomet- 
ric membership boundary. Thus, the present method 
of determining membership probabilities is likely to 
overestimate the membership probabilities for stars in 
the central region. 

In order to see the robustness of segregating the 
cluster member stars from the field stars, we created 
artificial clusters using the radial density profile and 
the mass function of NGC 1245 via the Monte-Carlo 
technique. The spatial distribution of each star was 
determined by randomly generated radial positions 
assuming the King profile whereas the positions in 
the CMD were determined from the luminosity func- 
tion assuming asymmetric Gaussian distribution for 
the colors at a given luminosity. We did not assume 
binary sequence explicitly but we took into account 
the binary effect by assuming three times larger a for 
the Gaussian distribution for the color distributions in 
the right side of the main sequence. We also generated 
the artificial field stars by assuming homogeneous spa- 
tial distribution and the CMD positions of the field 
stars which are located at P > 40' from the cluster 
centre. We derived the CMD position of each star by 
assigning magnitude and color randomly around each 
field star position in the CMD with AV = 0.125 and 
A.{B — V) = 0.025. We constrained the number of field 
stars by comparing the area where the field stars were 
selected and the total observed field of NGC 1245. 
In Fig. 5, we presented the spatial distributions and 
CMDs of an artificial cluster including artificial field 
stars. 

We ran 100 simulations to see that the 4-D mem- 
bership probability technique works correctly. Fig. 5 
shows the fractional frequency distributions of the 
cluster members and non-members as a function of 
the 4-D membership probability P. Most of the field 
stars have P less than ~ 0.2 with ^ 70% of field stars 
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Figure 5. Fractional frequency distribution of member- 
ship probability of the member and non-member stars. 
Solid and dotted lines represent member stars and non- 
member stars, respectively. 



have P = whereas ~ 50% of the cluster members 
have P larger than 0.2. As can be seen in Fig. 5, if 
we randomly select a star with P — 0.15, the prob- 
ability that this star is a field star is ~ 50%. But, 
if we randomly select a star with P = 0, the field 
star probability is '^ 70% whereas if we select a star 
with p > 0.4, that is less than ~ 10%. The reason for 
the relatively high number of member stars that have 
P = is the number of member stars is only 4% of the 
field stars. This makes P = for member stars which 
are located in the regions where the surface density 
is greatly dominated by field stars. A predominance 
of field stars makes Nc ~ and consequently P be- 
comes zero (see eqs (3) and (4)). Thus, the member 
stars that have P = are mostly located in the region 
where field stars are distributed with a wedge pattern 
in the CMD. 



4 STRUCTURE AND DYNAMICAL 
EVOLUTION OF OPEN CLUSTERS 

4.1 Radial Surface Density Profiles 

Similar to lChen. Chen and Shul (|200J), we defined an 
effective surface number density by summing up the 
membership probabilities. Pi, inside the observed area 

AS, 

f-— ^ f5) 

^- AS A ^^' 

where A is the correction factor for incomplete pho- 
tometry (Paper I). Because the position of maximum 
value of / {fmax) depends upon the size of AS, we 



first determined the position of fmax by calculating / 
for all the pixel positions {x,y) with AS = 1' x 1'. 
Then we determined the position of the cluster cen- 
tre by averaging the coordinates of the pixels where 
/ is larger than 0.9 fmax- We derived the radial sur- 
face number density f{R) by calculating / for concen- 
tric rings at a radial distance R with radial interval 
of dR = (l.l)"dPo where n represents the nth ring 
and dRo = l.'O is the radius of the central region. 
Although / is considered to be the surface number 
density corrected for field stars, there are some resid- 
uals, fbg due to field stars which have positive Pi owing 
to their CMD positions lying within the cluster mem- 
bership boundaries. We calculated ftg at P > 40' for 
NGC 1245 and and R > 35' for NGC 2506, respec- 
tively. Because Pi depends on the smoothing box and 
sampling resolution, we optimized them to minimize 
fbg. The adopted values of ftg are 0.20 ±0.62 for NGC 
1245 and 0.25 ±0.63 for NGC 2506. We su btracted fb g 
from / before fitting to the King profiles (|Kingjll962h . 
Fig. 6 and Fig. 7 show the logarithmic radial sur- 
face number density profiles of NGC 1245 and NGC 
2506 with the best-fit King profiles which can be ex- 
pressed as 



f{R)- 



fo 



1 ± {R/Rc 



(6) 



where fo is the central surface number density, Re is 
core radius. We derived fo and Re by least-squares fit- 
ting technique. The radius for the cluster extent was 
taken to be Rd , where the radial surface number den- 
sity profiles becomes constant. We derived Rd = 16' 
and 19' for NGC 1245 and NGC 2506, respectively 
(paper I). 

Since the effective radial surface number density 
near Rd is comparable to ftg, errors of ftg affect the 
derived fo and Re significantly. For NGC 1245, la 
variation of ftg gave rise to ~ 30% and ~ 50% varia- 
tions in fo and Re while they were reduced to ~ 10% 
to both fo and Re for NGC 2506. 

We examined the dependence of Re on the lu- 
minosity of stars because the spatial distribution of 
stars depends on the stellar mass if the clusters under 
consideration are dynamically relaxed. We plotted the 
effective radial surface number density profiles along 
with best-fitted King profiles to show this dependence 
in Fig. 6 and Fig. 7. As expected, the effective radial 
surface number density profiles of bright stars have 
smaller core radius than those of faint stars (see Table 
2). This is the very clear evidence of the mass segre- 
gation due to dynamical evolution of the clusters. 

One thing worth to mention is the departure of 
the observed surface number density from the King 
profile at P > 15', which is more pronounced in NGC 
1245. Although the difference between the observation 
and the King profile is within the Poisson errors, there 
is a possibility that they are caused by some real fea- 
tures such as tidal tails, especially in NGC 1245. We 
can see the morphology of these features in Fig. 9 
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Figure 6. King model fitting to the radial surface d ensity 
profile s for NGC 1245. Solid curve represents best fitlKinfJ 
111963) to the observed data. The error bars are Poisson 



Figure 7. King model fitting to the radial surface d ensity 
profile s for NGC 2506. Solid curve represents best fit lKing 
1119621 ) to the observed data. The error bars are Poisson 
errors. 



where we found that there are over-dense regions at 
R ~ 30' in the north-west of NGC 1245. However, the 
surface density of the over-dense region is about 3cr of 
the ftg , the morphology of the over-dense regions are 
much affected by the errors in ftg. 

4.2 Radial LFs 

We derived the luminosity function (LF) of the clus- 
ters using the membership probability defined in equa- 
tion 4. Because the number of stars at each magnitude 
interval, corrected for the field star contamination, is 
equal to the effective number of stars N^jf , it can be 
expressed as 



Table 2. King model fitting parameters for NGC 1245 and NGC 2506 



cluster 


Magnitude range 


/o 


Re 


NGC 1245 


y < 17 


8.3 


2.3 




17 < y < 20 


4.3 


4.5 




20 < y < 23 


3.1 


5.0 




y < 23 


14.7 


3.6 


NGC 2506 


y < 17 


14.3 


2.0 




17 < y < 20 


15.2 


3.7 




20 < y < 23 


10.4 


6.0 




y < 23 


33.5 


4.0 
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Figure 8. Radial dependence of the luminosity functions 
of NGC 1245 and NGC 2506. Clusters are divided into 
three regions, R < _Rc, Re < R < R\/2y ^nd ^1/2 < R < 
Rci , from the left to the right. The error bars are Poisson 
errors. 
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Figure 9. Log scale density contour map for NGC 1245. 
The inner contour represents /o/2 and the outer contour 
indicates 3cr level of background density fluctuation. Each 
arrows represent the Galactic coordinate (1, b) and the 
proper motion vector (PM). 



Neff{m) 



T.Pi{m) 

A(m) 



(7) 



where Pi{m) and A(m,) are the membership probabil- 
ity and incompleteness correction factor, respectively, 
for stars with a magnitude interval (m— 0.5 ~ m-f 0.5). 
Then, the LFs at a radial distance R from the cluster 
centre can be determined based on N^ffirn). 

We considered three regions, R < Re, R < R1/2, 
and Ri/2 < R < Rci, to derive Radial dependence of 
LFs. The half number radius R1/2 is the radius where 
the number of stars becomes half of the total number 
of stars within R^i- They are R1/2 ~ 7.'1 and 7.'8 for 
NGC 1245 and NGC 2506, respectively. 

Fig. 8 shows the radial LFs derived for the two 
open clusters, NGC 1245 and NGC 2506. The general 
trend of the LFs of the two open clusters under con- 
sideration is a gradual change of the peak positions; 
the LFs peak closes to the bright end of the inner re- 
gion {R < Re), whereas the peak closes to the faint 
end in the outer region {R1/2 < R < Rci)- This means 
that fainter stars are likely to be found in the outer 
regions of the clusters, i.e., there is a clear evidence of 
mass segregation in the two clusters, NGC 1245 and 
NGC 2506. 

Similar sign atures of mass segregation were re- 
ported fo£^Il]_jMathicu 1984), Pleiades and Prae- 
sepe (iRaboud fc Mermilliod 1998a .b.) NGC 2099 
(iKalirai et al. 112001b) and NGC 6819 (JKalirai et al. I 
l2001al : iKang fc AnnlboO^ )). The reason for the ob- 
served mass segregation seems to be the dynamical 



evolution that leads to the energy equilibrium because 
these open clusters are old enough to be rela:xed. 



4.3 Halo Structures 

To investigate the morphology and halo structure of 
the clusters, we plotted 2 dimensional effective surface 
number density distributions of NGC 1245 and NGC 
2506 in Fig. 9 and Fig. 10, respectively. We plotted two 
contour levels. One corresponds to the effective sur- 
face number density of /o/2 showing the representive 
morphology of the cluster and the other corresponds 
to the 3cr of the background density (fbg) showing the 
outer morphology of the cluster which can be detected 
above the background fluctuations. As can be seen in 
Fig. 9 and Fig. 10, NGC 1245 is more elongated than 
NGC 2506. There seems to be no relation between the 
direction of the major axis of NG C 1245 and the direc- 
tion of the cluster proper motion (jLoktin fc Beshenovl 
I2OO3I ) . We plotted the direction s of cluster proper mo- 
tion iLoktin fc Beshenovl l2003h along with the direc- 
tions of the Galactic latitude and longitude for com- 
parison. There are some features resembling tidal tails 
in the north-west of NGC 1245. The surface number 
densities of these features are comparable to the 3cr 
of the fbg which is ~ 1/8 times the central surface 
density (/o). 
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Figure 10. Log scale density contour map for NGC 2506. 
The inner contour represents /o/2 and the outer contour 
indicates Scr level of background density fluctuation. Each 
arrows represent the Galactic coordinate (1, b) and the 
proper motion vector (PM). 



5 SUMMARY AND CONCLUSIONS 

We proposed a new method for determining the clus- 
ter membership probability using the spatial positions 
together with the positions in the CMD. Owing to the 
ability to assign the membership probability for all the 
stars in the cluster field, we could trace the morphol- 
ogy of the clusters and derive the radial surface density 
profiles corrected for the field star contamination. 

The validity of the present method of 4-D 
membership probability is confirmed by compar- 
ing our membership probabilities with those from 
prope r motion studies in the central region of NGC 
2506 IjChiu fc va n Altcna 1981). However, our method 
shows some bias to higher probability for stars in the 
central regions of the cluster. This bias is caused by 
assigning high probabilities to stars in the central re- 
gions of the cluster if they are located within the 
photometric membership boundary in the CMDs. We 
also tested the robustness of the segregation of non- 
member stars from the member stars by the Monte- 
Carlo simulations and found that the 4-D method as- 
signs membership probability quite convincingly. The 
probability of a star with P > 0.4 is a real cluster 
member is about 90%. However, for better under- 
standing of the structure of open clusters, we need 
kinematic data. 

The two old open clusters NGC 1245 and NGC 
2506 show an imprint of the dynamical evolution in 
their LFs and radial surface density profiles. The LFs 
of the inner region of the clusters are dominated by the 
giant and upper main-sequence stars, whereas those of 



the outer region of the clusters are dominated by low- 
mass stars. This effect is more pronounced in NGC 
2506. Some fraction of low mass stars is expected to 
be evaporating from the cluster to create the tidal 
halo surrounding the cluster. The over-dense regions 
revealed in the effective number density contour map 
of NGC 1245 suggest the presence of the tidal halo at 
R > Rci in the north-west of NGC 1245. The detec- 
tion of the over-dense regions at _R > Rd was made 
possible due to our method of membership probabil- 
ity which took into account spatial positions together 
with the positions in the CMD, in addition to using 
wide and deep photometry of the two clusters (Paper 

I)- 

We found that the cluster core radius derived 
from bright stars is smaller than that from the faint 
stars. The dependence of core radius on the luminosity 
of stars is a strong evidence of the mass segragation 
in NGC 1245 and NGC 2506. Mass segregation and 
evaporation of the low mass stars are natural conse- 
quences of the dynamical evolution in relaxed clusters. 
Stars exchange their kinetic energy through encoun- 
ters. Owing to the equipartition of the kinetic energy, 
high-mass stars have relatively low velocity, whereas 
low mass stars may have velocity that is sufficient to 
escape from the cluster. If the evaporated stars are 
tidally disturbed by the Galactic disc or giant molec- 
ular clouds during their orbital motion, they are likely 
to be stretched along the orbit (|Chumak et al.ll20ld ) 
to create tidal tails. A signature of tidal tail seems to 
be observed in NGC 1245 but the signal is too low to 
be confirmed. 
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